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Abstract 



The left-invariant sub-Riemannian problem on the group of motions 
^^ . of a plane is considered. Sub-Riemannian geodesies are parametrized by 

Jacobi's functions. Discrete symmetries of the problem generated by re- 
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1 Introduction 

Problems of sub-Riemannian geometry have been actively studied by geometric 
control methods. One of the central and hard questions in this domain is a 
description of cut and conjugate loci. Detailed results on the local structure of 
conjugate and cut loci were obtained in the 3-dimcnsional contact case [1,4]. 
Global results are restricted to symmetric low-dimensional cases, primarily for 
left-invariant problems on Lie groups (the Heisenberg group [6,20], the growth 
vector {n,n{n + l)/2) [8-10], the groups S0(3), SU(2), SL(2) and the Lens 
Spaces [5]). 

The paper continues this direction of research: we start to study the left- 
invariant sub-Riemannian problem on the group of motions of a plane SE(2). 
This problem has important applications in robotics [7] and vision [11]. On the 
other hand, this is the simplest sub-Riemannian problem where the conjugate 
and cut loci differ one from another in the neighborhood of the initial point. 

The main result of the work is an upper bound on the cut time <cut given 
in Theorem l6.4l we show that for any sub-Riemannian geodesic on SE(2) there 
holds the estimate tcut < t, where t is a certain function defined on the cotangent 
space at the identity. In a forthcoming paper [19] we prove that in fact icut = t. 
The bound on the cut time is obtained via the study of discrete symmetries of 
the problem and the corresponding Maxwell points — points where two distinct 
sub-Riemannian geodesies of the same length intersect one another. 

This work has the following structure. In Section [5] we state the problem 
and discuss existence of solutions. In Section [3] we apply Pontryagin Maxi- 
mum Principle to the problem. The Hamiltonian system for normal extremals 
is triangular, and the vertical subsystem is the equation of mathematical pen- 
dulum. In Section 0] we endow the cotangent space at the identity with special 
elliptic coordinates induced by the flow of the pendulum, and integrate the nor- 
mal Hamiltonian system in these coordinates. Sub-Riemannian geodesies are 
parametrized by Jacobi's functions. In Section [5] we construct a discrete group 
of symmetries of the exponential mapping by continuation of reflections in the 
phase cylinder of the pendulum. In the main Section[6]we obtain an explicit de- 
scription of Maxwell strata corresponding to the group of discrete symmetries, 
and prove the upper bound on cut time. This approach was already success- 
fully applied to the analysis of several invariant optimal control problems on Lie 
groups [13-18]. 



2 Problem statement 

The group of orient 
resented as follows: 

SE(2) = i sin6l cos6' y \ \ 6 e S^ =R/{2ttZ), x,y e 



The group of orientation-preserving motions of a two-dimensional plane is rep- 
resented as follows: 



cosy 


— smy 


X 


sin 6 


cos 6* 


y 
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The Lie algebra of this Lie group is 

e(2) = span(^2i - -E12, S13, £^23), 

where Eij is the 3x3 matrix with the only identity entry in the i-th row and 
j'-th column, and all other zero entries. 

Consider a rank 2 nonintegrable left-invariant sub-Riemannian structure on 
SE(2), i.e., a rank 2 nonintegrable left-invariant distribution A on SE(2) with a 
left-invariant inner product (•, •) on A. One can easily show that such a structure 
is unique, up to a scalar factor in the inner product. We choose the following 
model for such a sub-Ricmannian structure: 

Ag = span(^i((7), 6 (<?)), (6:, 6) = '^U' h 3 = 1,2, 

6(9) =9^13, S,2{q)=q{E2i- E12), 
and study the corresponding optimal control problem: 

q^uiii{q)+U2i2{q), qeM = SE{2), u ^ {ui,U2) e R^ , (2.1) 

g(0) = <?o = ld, q{ti)^qi, 

£ — \/ul + ul dt -^ min . 

In the coordinates (x, y, 9), the basis vector fields read as 

a = COS 0— -I- sine— , 6 = TT7T, (2.2) 

ox ay oO 

and the problem takes the following form: 

x — uiC0s9^ y — uisa\9, 9 — U2, (2-3) 

q=ix,y,9)eM = RlyXSl, u = (ui, U2) e M^, (2.4) 

g(0) = 90 = (0,0,0), q{h)^q, = ixi,yi,9i), (2.5) 

i^ Jul + uldt-^ min . (2.6) 

Admissible controls u(-) are measurable bounded, and admissible trajectories 
g(-) are Lipschitzian. 

The problem can be reformulated in robotics terms as follows. Consider 
a mobile robot in the plane that can move forward and backward, and rotate 
around itself (Reeds-Shepp car). The state of the robot is described by coordi- 
nates (x, y) of its center of mass and by angle of orientation 9. Given an initial 
and a terminal state of the car, one should find the shortest path from the initial 
state to the terminal one, when the length of the path is measured in the space 
{x,y,9), see Fig.[Tl 

Cauchy-Schwarz inequality implies that the minimization problem for the 
sub-Riemannian length functional (j2.6p is equivalent to the minimization prob- 
lem for the energy functional 

1 Z"*^ 

J = - / {ul + ul) dt -> min (2.7) 

2 Jo 




Figure 1: Problem statement 



with fixed ti. 

System ()2.1|) has full rank: 



d d 

6 = Ki,6]=sin0— -cos^— , 

span(a(g),e2(g),6('?)) = TgM V q G M, 



(2.8) 
(2.9) 



so it is completely controllable on M. 

Another standard reasoning proves existence of solutions to optimal control 
problem (|2.1[) . (|2.5p . (|2.7[) . First the problem is equivalently reduced to the time- 
optimal problem with dynamics (|2.ip . boundary conditions (j2.5p . restrictions on 
control uf + U2 < 1, and the cost functional ti — > min. Then the state space of 
the problem is embedded into M.^ (see e.g. [17]), and finally Filippov's theorem [3] 
implies existence of optimal controls. 



3 Pontryagin Maximum Principle 

We apply the version of PMP adapted to left-invariant optimal control problems 
and use the basic notions of the Hamiltonian formalism as described in [3]. 
In particular, we denote by T*M the cotangent bundle of a manifold M, by 
TT : T*M -^ M the canonical projection, and by /i S Vec(T*M) the Hamiltonian 
vector field corresponding to a Hamiltonian function h S C°"{T*M). 

Consider linear on fibers Hamiltonians corresponding to the vector fields S,i '■ 

/i,(A) = (A,6(g)), g = 7r(A), XeT*M, i = 1,2,3, 

and the control-dependent Hamiltonian of PMP 

Ki\)='^{ul + ul)+uihi{X)+U2h2{X), XeT*M, weM^ ly e {-1,0}. 



Then the Pontryagin Maximum Principle [3, 12] for the problem under con- 
sideration reads as follows. 

Theorem 3.1. Let u{t) and q{t), t G [0,ti], be an optimal control and the 
corresponding optimal trajectory in problem \2.1^ . ^2.5^ . \2.'7^ . Then there exist 
a Lipschitzian curve At € T* M , 7r(At) — q{t), t G [0,ti], and a number v G 
{—1, 0} ]or which the following conditions hold for almost all t G [0, ti]: 



At = /i;:(t)(At) = ui{t)hi{Xt) + u2{t)h2{Xt), 


(3.1) 


Kit)(M) =max/i5:(At), 


(3.2) 


(^At)^O. 


(3.3) 



Relations (|2.8p , (|2.9p mean that the sub-Riemannian problem under consid- 
eration is contact, thus in the abnormal case 1^ = the optimal trajectories are 
constant. 

Consider now the normal case v = —1. Then the maximality condition p.2p 
implies that normal extremals satisfy the equalities 

Ui{t) = ht{Xt), i = 1,2, 

thus they are trajectories of the normal Hamiltonian system 

A = i?(A), A G T*M, (3.4) 

with the maximized Hamiltonian H = {hf + h2)/2. 
In view of the multiplication table 

[a,6]-6, [a,6] = o, [6,6] = a, 

system (|3.4p reads in coordinates as follows: 

hi = -h2h3, h2 = hihs, ha = /i2^3, (3.5) 

i = hiCos6, y = hism9, 6 = h2. (3.6) 

Along all normal extremals we have H = C > 0; moreover, for non-constant 
normal extremal trajectories C > 0. Since the normal Hamiltonian system (|3.5p . 
(|3.6p is homogeneous w.r.t. (/ii,/i2), we can consider its trajectories only on 
the level surface H = 1/2 (this corresponds to the arc-length parametrization 
of extremal trajectories), and set the terminal time ti free. Then the initial 
covector A for normal extremals At = e*^(A) belongs to the initial cylinder 

c^T;Mn{H{\) = i/2}. 

Introduce the polar coordinates 

hi — cos a, /i2 = sin a, 



then the initial cyhnder decomposes as C = 5*^ x R/j.,, where S^ = R/(27rZ). In 
these coordinates the vertical part (|3.5p reads as 



a — h^, h3 — -sm2a, {a,h3)£C. (3.7) 

In the coordinates 

7 = 2a + 7re 25"^ =M/(47rZ), c = 2/i3 e M, 

system (|3.7p takes the form of the standard pendulum 

Af^c, c = - sin 7, (7, c) e C ^ (25^ ) x M^. (3.8) 

Here 25*1 = R/(47rZ) is the double covering of the standard circle 5*^ = R/(27rZ). 
Then the horizontal part (|3.6p of the normal Hamiltonian system reads as 



i; = sin — COS0, y = sin — sin 6*, 9 — —cos—. (3.9) 

Summing up, all nonconstant arc-length parametrized optimal trajectories in 
the sub-Riemannian problem on the Lie group SE(2) are projections of solutions 
to the normal Hamiltonian system ()3.8|) . (13. 9|) . 



4 Exponential mapping 

The family of arc-length parametrized normal extremal trajectories is described 
by the exponential mapping 

Exp : N -> M, iV = C X R+, 

Exp(zy) = Exp(A, i) = TT o e}^{X) = n{Xt) = q{t), 
,^ = iX,t) = {-f,c,t)eN. 

In this section we derive explicit formulas for the exponential mapping in special 
elliptic coordinates in C induced by the flow of the pendulum p.8p . The general 
construction of elliptic coordinates was developed in [13, 14, 17], here they are 
adapted to the problem under consideration. 

4.1 Decomposition of the cylinder C 

The equation of pendulum (j3.8|) has the energy integral 

c2 

^ = — -cos7e [-l,+oo). (4.1) 



Consider the following decomposition of the cylinder C into disjoint invariant 
sets of the pendulum: 



i=l 

Ci = {A e C 
C2 = {A e C 
C3 = {A e C 
C4 = {A e C 
C5 = {A e C 



(4.2) 



i?e(-i,l)}, 
Ee (l,+oo)}, 

E=l, c^Q], 

E = -l} = {(7, c) e C I 7 = 27rn, c = 0}, 

E = l, c = 0} = {(7, c) e C I 7 = TT + 27r7i, c = 0}. 



0,1, 



Here and below we denote by n a natural number. 
Denote the connected components of the sets C^: 

Ci - uLo^i, C\ = {(7,c) e Ci I sgn(cos(7/2)) = (-1)*}, z 

C2 = C+ U C2- , C± = {(7, c) e C2 I sgn c = ±1}, 

C3 = uLo(Q+UCr), 

C^± = {(7,c)eC3|sgn(cos(7/2)) = (-l)\ sgnc-il}, * 



0,1, 



C5 = Uj^qC's, 



Q = {(7, c) e C I 7 = 2m, c = 0}, i = 0, 1, 



C^ 



{(7, c) e C I 7 = TT + 27ri, c = 0}, i == 0, 1. 



Decomposition (14. 2p of the cylinder C is shown at Fig. [21 
c 
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Figure 2: Decomposition of the cylinder C 



4.2 Elliptic coordinates on the cylinder C 

According to the general construction developed in [17], we introduce elliptic 
coordinates {(f, k) on the domain Ci U C2 U C3 of the cylinder C, where A; is a 



reparametrized energy, and ^ is the time of motion of the pendulum (|3.8p . We 
use Jacobi's functions am((/j, fc), cn((p, fc), sn((^, A:), dn((/3, A:), E((y9, fc); moreover, 
i^(fc) is the complete elliptic integral of the first kind [21]. 
If A = (7, c) e Ci, then: 



2 

If A = (7,c) e C2, then 



sin — = 5iA;sn((^, /c), Si = sgn 003(7/2), 

cos- = sidn(v3. A:), 
^=fccn(^,/c), (^e[0,47^(fc)]. 



7 
sin- = S2Sn((^/fc,fc), S2=sgnc, 

7 
cos — = cn((^/fc, /c), 

I = (s2/fc) dn(^/fc. A:), v? € [0, 4A;X(A;)]. 

If A = (7, c) G C3, then: 

A: = 1, 

7 
sin — = S1S2 tanh(p, si = sgncos(7/2), S2 — sgnc 

7 
cos — — Si/ cosh (y9, 

- = 82/ cosh 93, (^ G (— cx3,+oo). 

4.3 Parametrization of extremal trajectories 

In the elliptic coordinates the flow of the pendulum ()3.8|) rectifies: 

(/3 = 1, fc = 0, A= (93,A;) euf^iCi, 



this is verified directly using the formulas of Subscc. 14.21 Thus the vertical sub- 
system of the normal Hamiltonian system of PMP p.Sp is trivially integrated: 
one should just substitute Lpt = Lp + t, k = const to the formulas of elliptic co- 
ordinates of Subsec. 14.21 Integrating the horizontal subsystem p. 91) . we obtain 
the following parametrization of extremal trajectories. 



If A = {(f, k) E Ci, then ipt — ^ + t and: 

cos Of — cuif en ipt + snip sn ipt , 

sin^f — si{s'nipcTiipt — cinp sn ipt) , 

Of — si{£Lniip — amipt) (mod 27r), 

Xt = {si/k)[cnip{dnip - dmpt) +snip{t + £(93) - E(i^t))], 

yt = {l/k)[snip{dnip - dnipt) - cnip(t + E{(p) - E((/Jt))]. 

In the domain C2, it wiU be convenient to use the coordinate 

%l)^tp/k, iPf = ipt/k ^ i; + t/k. 

If A e C2, then: 

cos Ot — k sn -0 sn if^t + dml) dn ipt , 

sinOt = fc(sn'0 dnipt — dmp snipt), 

Xt = S2fc[dnV'(cnV' - cnV't) + snip{t/k + E(V') - E(V't))], 

yt = S2[fc^ sn-0(en-0 — en-0t) — dmp{t/k + E{ip) — Eiipt))]- 

If Ae C3, them 

cos 9t = l/(cosh(^coshiy9t) + tanh (^9 tanh (^94 , 

sin^f — Si(tanhiy9/ eoshi^s^ — tanhipt/ coshes), 

Xt = siS2[(l/ cosh 1^9) (1/ eosh(y9 — 1/ eoshi^Si) + tanhiy9(i + tanh(^ — tanhipt)], 

yt = S2[tanh 1^9(1/ eosh(y9 — 1/ cosh(y9f) — (1/ cosh 1^9) (i + tanh(y9 — tanhiy94)]. 

In the degenerate cases, the normal Hamihonian system p.Sp . (|3.9|1 is easily 
integrated. 

If A e C4, then: 

Ot = -Sit, Xt =0, yt ^ 0. 

If Ae C5, them 

9t = 0, Xt = t sgnsin(7/2), yt = 0. 



It is easy to compute from the Hamiltonian system p.Sp . p.9p that pro- 
jections (xt,yt) of extremal trajectories have curvature k = — cot(7t/2). Thus 
they have inflection points when cos(7t/2) = 0, and cusps when sin(7i/2) = 0. 
Each curve {xt,yt) for A G U^^^Ci has cusps. In the case A G Ci U C3 these 
curves have no inflection points, and in the case A € C2 each such curve has 
inflection points. Plots of the curves (a:^, yt) in the cases A € Ci U C2 U C3 are 
given respectively at Figs. [31 IH [5l 

In the cases A G C4 and A G C5 the extremal trajectories qt are respectively 
Riemannian geodesies in the circle {x = y = 0} and in the plane {9 = 0}. 



10 





Figure 3: Non-inflexional trajec- 
tory: A e Ci 



Figure 4: Inflexional trajectory: A € 
C2 




Figure 5: Critical trajectory: A G C3 
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5 Discrete symmetries and Maxwell strata 

In this section we continue reflections in the state cyUnder of the standard 
pendulum to discrete symmetries of the exponential mapping. 

5.1 Symmetries of the vertical part of Hamiltonian system 
5.1.1 Reflections in the state cylinder of pendulum 

The phase portrait of pendulum ([3 



m (|,3.8| 


admits the following reflections: 


(7,c)- 


* (7,-c), 


(7,c)- 


^ (-7,c), 


(7, c) - 


^ (-7,-c), 


(7, c) - 


>(7 + 27r,c), 


(7, c) - 


>(7 + 27r,-c), 


(7, c) - 


>(-7 + 27r,c), 


(7,c)- 


>(-7 + 27r,-c). 



These reflections generate the group of symmetries of a parallelepiped G = 
{Id, £^, . . . je"^}. The reflections e'^, e'^, e^ preserve direction of time on trajec- 
tories of pendulum, while the reflections e^, e^, e^, e^ reverse the direction of 
time. 



5.1.2 Reflections of trajectories of pendulum 

Proposition 5.1. The following mappings transform trajectories of pendu- 
lum l\3.8\ to trajectories: 



e* : 5 = {(7„c.)|,se[0,i]}K-<5^ = {(7:,c:)|sG[0,i]}, z = l,...,7, (5.1) 
where 

(7s,c^) = (-7t_s,ct-.), 
(7',c^) = (-7„-c3), 
(7f,4)-(7.+27r,c,), 
(7f,c^) = (7t_, + 2^,-Ct_,), 
(7f,c^) = (-7i_,+27r,Q_,), 
(7j,cI) = (-7. + 2^,-c,). 



Proof. The statement is verified by substitution to system (|3.8p and differenti- 
ation, n 



The action (|5.ip of reflections e* on trajectories 5 of the pendulum p.Sp is 
illustrated at Fig. [H 
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Figure 6: Reflections e* 



(5* of trajectories of pendulum 



5.2 Symmetries of Hamiltonian system 

5.2.1 Reflections of extremals 

We define action of the group G on the normal extremals As = e''^(Ao) G T*M, 
s S [0,t], i.e., solutions to the normal Hamiltonian system 



as follows: 



7s 



Cs, 



■sm7s 



sin y Xi {qs ) - cos y ^2 (Qs ) 



(5.2) 
(5.3) 



e' : {As|sG[0,f]}^{A:|se[0,t]}, 



,7, 



(5.4) 
(5.5) 



Here \l is a solution to the Hamiltonian system (|5.2p . (|5.3p . and the action of 
reflections on the vertical coordinates (7s, Cg) was defined in Subsec. l5.ll The 
action of reflections on the horizontal coordinates {xs,ys,ds) is described as 
follows. 

Proposition 5.2. Let qs — {xs,ys,Os), s G [0,i\, be a normal extremal trajec- 
tory, and let ql — {x\,y].,9\), s € [0,t], be its image under the action of the 
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reflection e' as defined by ^5.4h ^KM- Then the following equalities hold: 

(1) &i — 9t — Ot_s, 

x] = cos0tixt - xt-s) + sin9t{yt - Vt^s), 
yl = sin 6t{xt - Xts) - cos 9t{yt - yt-s), 

(2) 02 = dt — Ot-s, 
Es = -cos9t{xt - xt-s) -sm9t{yt - yt-s), 



(3) 



vl^- sin 
^1= Os, 



9t{xt - xt-s) + cos6t{yt ~ yt~s), 



Xg Xs, 

(4) ei = -Os, 

4 _ _ 
Xg — Xg, 

yt = ys, 

(5) ei = et.s-Ot, 



x\ 

y' 



(6) 



y\ 
(7) 0? = -Qs, 



5 — "J — s "ti 

■\ = cos0t{xt-s - xt) + siTLet{yt-s -yt), 
■f = - sm9t{xt-s - Xt) + cos9t{yt-s - yt), 

6 = (^t-s - 9t, 

xf = cos6'f(xf -Xt-s) + s[n9t{yt - yt-s), 
yt = -sm9t{xt -xt_s) +cos9t{yt -yt-s), 



X 



yl = -ys- 



Proof. We prove only the formulas for 9l and xl since all other equalities are 
proved similarly. 

By Proposition l5.ll we have 7^ = jt-s- Then we obtain from 



and 



ps 1 pS pt — S 

91 = — COS —^dr = — / cos dr = cos — dp = 9t — 9t-s 

'Jo 2 Jo 2 Jt 2 ^ * * ^ 



x] — sin — ^ cos 9l dr = / sin cos{9t — 9t-r) dr 

Jo 2 Jo 2 

1*1 — S pt — S 

~ — COS 9t / sin -^ cos 9p dp — sin 9t / sin -;^ sin 9p dp 
Jt ^ Jt ^ 



cos 6t{xt - Xt-s) + sin 9t{yt -yt-s)- 
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D 



The action of reflections e* on curves {xs,ys) has a simple visual meaning. 
Up to rotations of the plane {x,y), the mappings e^ , e^, s^ are respectively 
reflections of the curves {{Xs,ys) \ s e [0,t]} in the center of the segment I 
connecting the endpoints (xo,j/o) and {xt,yt), in the middle perpendicular to 
I, and in I itself (see [14,17]). The mapping e^ is the reflection in the axis y 
perpendicular to the initial velocity vector (cos 0o 7 sin 6*0). The rest mappings 
are represented as follows: e'"^"' = e^ o e'\ i = 1,2, 3. 



5.2.2 Reflections of endpoints of extremal trajectories 

We define action of reflections in the state space M as the action on endpoints 
of extremal trajectories 



e' : M ^ M, 



qt ^-^ Qt, 



(5.6) 



see (|5.4p . (|5.5I) . By virtue of Propos. 15. 2[ the point ql depends only on the 
endpoint qt, not on the whole trajectory {qs \ s G [0,t]}. 

Proposition 5.3. Let q = {x,y,0) E M, q' = e'{q) = {x\y\e') G M. Then: 



X cos + y sm 0, x sm — y cos 0, 0) 

—x cos 9 -- y sin 9, —x sin 9 + y cos 6 

~x,-y,9), 

~x,y,~9), 

—x cos 9 — y sin 9, x sin 9 — y cos 6, - 

X cos 9 + y sin 9, —x sin + y cos 6, - 

x,-y,-9). 



{x' 


y' 


e') = { 


ix' 


y' 


e')^i 


{x' 


y' 


03) = ( 


{x' 


y' 


9') = { 


{x^ 


y' 


^') = ( 


{x^ 


y' 


^') = ( 


{x' 


y' 


9') = { 



Proof. It suffices to substitute s = to the formulas of Proposition [5?2] 



D 



5.3 Reflections as symmetries of exponential mapping 

Define action of the reflections in the preimage of the exponential mapping: 

e^:N^N, e' : ,. ^ ij,c,t) ^ ,.^ ^ {Y,c\t), (5.7) 

where (7, c) — (70, co) and (7*,c*) — {jo,Cq) are the initial points of the corre- 
sponding trajectories of pendulum (7^,0^) and (7^,0^). The explicit formulas 
for (7*, c*) are given by the following statement. 

Proposition 5.4. Let v = (A,<) = (7, c, i) e N , v^ ^ e'{iy) = {X\t) = 
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{-/\c\t) e N. Then: 

(7\c^) = {lt,-ct), 
(7^c^) = (-7t,ct), 
(7^c3) = (-7,-c), 
(7^c4) = (7 + 2vr,c), 

(7^c5) = (7t + 2^,-Q), 
(7^c6) = (-7, + 27r,Q), 
(7',c^) = (-7,-c). 
Proof. Apply Proposition 15. II with s = 0. D 

Formulas (|5.6p . (|5.7p define the action of reflections e* in the image and 
preimage of the exponential mapping. Since the both actions of e' in M and N 
are induced by the action of e' on extremals As ()5.4p . we obtain the following 
statement. 

Proposition 5.5. For any i = 1, . . . ,7, the reflection e* is a symmetry of the 
exponential mapping, i.e., the following diagram is commutative: 

Exp Exp 

N >■ M v\ *-9 





6 Maxwell strata corresponding to reflections 

6.1 Maxwell points and optimality 
of extremal trajectories 

A point qt of a sub-Riemannian geodesic is called a Maxwell point if there exists 
another extremal trajectory qs ^ <ls such that qt — Qt for the instant of time 
i > 0. It is well known that after a Maxwell point a sub-Riemannian geodesic 
cannot be optimal (provided the problem is analytic). 

In this section we compute Maxwell points corresponding to reflections. For 
any z = 1, . . . , 7, define the Maxwell stratum in the preimage of the exponential 
mapping corresponding to the reflection e* as follows: 

MAX* ^{v={\t)(^N\\^\\ Exp(A, t) = Exp(A\ i)}. (6.1) 

We denote the corresponding Maxwell stratum in the image of the exponential 
mapping as 

Max' = Exp(MAX*) C M. 

li V ^ {\,t) £ MAX*, then qt = Exp(i') e Max* is a Maxwell point along the 
trajectory qs — Exp(A, s). Here we use the fact that if A ^^ A*, then Exp(A, s) ^ 
Exp(A*,s). 
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6.2 Multiple points of exponential mapping 

In this subsection we study solutions to the equation q = q^ ^ where g* — e*(g), 
that appears in definition (|6.ip of Maxwell strata MAX*. 

The following functions are defined on M = R^^y x Sg up to sign: 

^ . S r. 9.9 

Ri = V cos X sm — , Ro ~ X cos — \- x sm — , 

i w 2 2 2 2 

although their zero sets {Ri — 0} are well-defined. In the polar coordinates 

a; = pcosx, y = psinx, 

these functions read as 

7?i =psm I X- 2 j ' -R2=pcoslx-2 

Proposition 6.1. (1) q^ = q <^ Ri{q) = 0. 

(2) q^=q ^ R2{q)=0- 

(3) q^ = q <^ x = y ^0. 

(4) q"^ =g ^ sine = a; = 0. 
{5)q^^q 4^ e = n or {x,y, 9) = (0,0,0). 

(6) q^ ^q ^ 9 = 0. 

(7) q"^ =q ^ sine = y = 0. 

Proof. We prove only item (1), all the rest items are considered similarly. By 
virtue of Proposition l5.3[ we have 

1 I X cos 9 + y sin 9 = x In sin ^ sin f Y — I) = 

I a; sin e — y cos 9 = y \p cos | sin (x ~ §) =0 

4^ psin(x-7r) <^ Ri{q) = 0. 



n 

Proposition 16.11 implies that all Maxwell strata corresponding to reflections 
satisfy the inclusion 

Max^ C{qe M \ Ri{q)R2{q) sm9 = 0}. 

The equations Ri{q) = 0, i = 1,2, define two Moebius strips, while the equation 
sine = determines two discs in the state space M = M^ x Sg, see Fig. [T] 

By virtue of Propos. 16. 1[ the Maxwell strata Max^, Max , Max are one- 
dimensional and are contained in the two-dimensional strata Max"'^, Max^, Max^, 
Max . Thus in the sequel we restrict ourselves only by the 2-dimensional strata. 
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Figure 7: Surfaces containing Maxwell strata Max* 
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6.3 Fixed points of reflections 

in preimage of exponential mapping 

In this subsection we describe solutions to the equations A = A^ essential for 
explicit characterization of the Maxwell strata MAX', see (I6.ip . 

From now on we will widely use the following variables in the sets Ni, i = 
1,2,3: 

iy^iX,t)eNi ^ T^iip + ipt)/2, P^t/2, 

v = {\t)eN2 => T=(^ + ^t)/(2fc), p = t/{2k), 

iy=iX,t)eN3 ^ T={ip + ^t)/2, P = t/2. 

Proposition 6.2. Let (A,i) € N, £*(A,t) = (A\ t) G N. Then: 

(1)A^ = A ^ 1"^^=°' '^^^' ^ 

I is impossible for A G C2 U C3, 

(2)A^ = A ^ 1^"^='^' AeCrUC,, 
\t = AeCs, 

(3) A^ = A is impossible, 

(4) A^ = A <^ 



g ^ J is im,possible for A G Ci U C3, 



cnr =0, A e C2. 

Proof. We prove only item (1), all other items are proved similarly. 
By Propos.[53 if A S Cj, then A^ G Cj, i = 0, 1. Moreover, 

, 1 , \lt — 1 \snwt ~ snw 

\^ = \ ^ {' ^ -^ { ^ ^ <:^ cnT=0. 

I —Ct — c I — en (y9f = en iy9 

If A G C2 , then A^ G C^ , thus the equality A^ = A is impossible. 
Similarly, if A G C^-^ , then A^ G C^ , i = 0, 1, and the equality A^ = A is 
impossible. D 

6.4 General description of Maxwell strata 
generated by reflections 

We summarize our computations of the previous subsections. 
Theorem 6.1. Let v = (A,t) G U^^^iVj and qt = (xt,yt,Ot) = Exp(iy). 

(l).GMAX^ ^ f^i('Z*) = 0, cnr^O, /or A G C„ 

\ii'i((Zt)-0, /or A GC2UC3. 

(2)i/gMAX2 ^ i?2(gt)==0, suTT^O. 
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(3)i^GMAX5 ^ 0t=Tr or {xt,yt,Ot)^ (0,0,0). 
|6't = 0, cnr 7^0 /or Ae Cj. 



Proof. Apply Propositions 16.11 and 16.21 D 

6.5 Complete description of Maxwell strata 

We obtain bounds for roots of the equations Ri{qt) — 0, sin 6*4 = that appear 
in the description of MaxweU strata given in Th. 16.11 

We use the foUowing representations of functions along extremal trajectories 
obtained by direct computation. 

If AeCi, then 

sin^t == -si •2cnpi3dnT/A, (6.2) 

cos{et/2) = S3 ■ cnp/VA, (6.3) 

sin(6lt/2) = S4-6dnT/%/A, (6.4) 

i?i('Zt) = -S3-2(p-E(p))cnr/(fcVA), (6.5) 

i?2(gt) = -s4•2/2(p,fc)snT/(fc^/A), (6.6) 

f2{p,k) = k cnpQ — dnp{p — E{p)), 
A = l- k'^sn^p dn^T, 
S3 = ±1, S4 = ±1, si = -5354- 

If Ae C2, then 

sin Ot = -2fcfi dnp en r / A, (6.7) 

cos(6't/2) = S3 • dnp/VA, (6.8) 

sin(6't/2) = S4-fci3cnT/VA, (6.9) 

Riiqt) = S2S4 • 2{p - E(p)) dn r /VA, (6.10) 

R2{qt) = S2S4 • 2kfi{p, k) snr /%/A, (6.11) 

h{p,k) = cnp(E(p)-p)- dnpfi, (6.12) 

S3 = -S4 = ±1. 

Proposition 6.3. Let t > 0. 

(l) If X e Ci, then 9t = ^ p = 2Kn. 

{2)If\eC2,thenet^0 ^ (p = 2Kn or cut ^ 0). 

(3) If X (z C3, then 6t — is impossible. 

Proof. Apply (|6.4[) in item (1), (|6.9[) in item (2), and pass to the limit fc ^ 1 — 
in item (3). D 
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Proposition 6.4. Let t > 0. 

(1) IfXeCi, then 9t^n ^ p = K + 2Kn. 

(2) //A G C2, then 9t — n is impossible. 

(3) If X (z C3, then 6t — n is impossible. 

Proof. Apply (|6.3p in item (1), (|6.8p in item (2), and pass to the limit k ^ 1 — 
in item (3). D 

Lemma 6.1. For any k G (0, 1) and any p > we have p -~ E(p) > 0. 

rp rP 

Proof p - E(p) ^p- I dT?tdt^k^ I ST?tdt>0. D 

Jo Jo 

Proposition 6.5. Let t > 0. 

{l)IfXeCi,thenRi{qt)^0 ^ cnr=0. 

(2) // A G C2, then Ri{qt) — is impossible. 

(3) If X (z C3, then Ri{qt) = is impossible. 

Proof. Apply (|6.5p and Lemma [6.11 in item (1); (|6.10p and Lemma [6.11 in item 
(2); and pass to the limit fc ^ 1 — in item (3). D 

Lemma 6.2. For any k G (0, 1) and p > we have f2{p, k) > 0. 

Proof. The function /2 {p) has the same zeros as the function 52 (p) — /2 (p)/ dnp . 
But 52 (p) > forp > since 32(0) = and 33 (p) = fc^cn^p/dn^p > 0. D 

Lemma 6.3. For any k G [0, 1), the function fi{p) has a countable number of 
roots 

p^Pi{k), 71 G Z, 

p?=0, pr"(fc)--K(fc)- (6.13) 

The positive roots admit the bound 

p1{k)e{-K + 2Kn,2Kn), n G N, A:G(0, 1), (6.14) 

^5^(0) = 27m, n G N. (6.15) 

All the functions k 1— > Pi{k), n ^ "L, are smooth at the segment k G [0, 1). 

Proof. The function /i (p) has the same roots as the function gi (p) = /i (p)/ cnp . 
We have 

g[ip) = -dn^p/cn^p, (6.16) 

so the function gi{p) decreases at the intervals p G {—K + 2Kn, K + 2Kn), 
n G Z. In view of the limits 

51 (p) ~^ ioo as p -^ K + 2Kn ± 0, 
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the function gi (p) has a unique root p = pi at each interval p € {—K+2Kn, K+ 
2Kn), neZ. 

For p — 2Kn, n £ N, we have gi{p) = E(p) — p < 0, thus the bound (|6.14p 
follows. 

Further, equality (|6.15p follows since fi{p,0) = — sinp. 

Equalities (|6.13p follow since the function fi{p) is odd. 

By implicit function theorem, the roots Pi{k) of the equation gi{p) — are 
smooth in k since g'i{p) < when cnp ^ 0, see (|6.16p . D 

Corollary 6.1. (1) The first positive root of the function fi{p) admits the 
bound 

pl{k) e {K{k), 2K{k)), k e (0, 1). (6.17) 



(2) Ifpe{Q,pl), thenf,{p)<0. 

(3) lim Pi{k) = vr, lim Pi{k) = +cx). 

k — >+0 k — *1 — 

Plots of the functions K{k), pl{k), 2K{k) are given at Fig. [51 



Figure 8: Plots of the functions K{k) < 
pl{k) < 2K{k) 



Proposition 6.6. Let t > 0. 

{l)IfXeCi,thenR2{qt)^0 ^ snr=0. 

(2) // A e C2, then i?2(gt) = ^ (p = K(^) or snr = 0). 

{3)IfX£C3,thenR2{qt) = ^ t = 0. 



Proof. Apply (j6.6p and Lemma 16.21 in item (1), (j6.1ip and Lemma 
(2), and pass to the limit fc ^ 1 — in item (3). 

Lemma 6.4. If v e N1UN2UN3, then {xt,Vt,9t) ^ (0,0,0). 



in item 
D 
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Proof. The equality {xt,yt,Ot) = (0,0,0) is equivalent to {Ri{qt),R2{qt),0t) = 
(0,0,0). 

li 1/ e A^i, then the equalities Ri{qt) — 0, i?2(<Zt) = are equivalent to 
CUT =0, snr =0 fPropos. 16.51 16. 6p . which are incompatible. 

li v ^ N2U Ns, then the equality Ri (qt) =0 is impossible (Propos. 16. 5p . D 

On the basis of results of the previous subsections we derive the following 
characterization of the Maxwell strata. 

Theorem 6.2. (1) MAX^ (iNi = MAX^ nN2 = MAX^ (iNs = 0. 

(2) MAX^ niVi = MAX^ niVa = 0, 

MAX^ niVs = {i^eN2\p = p'l{k), snr 7^ 0}. 

(3) MAX^ niVi ^{vC,Ni\p = K + 2Kn}, 
MAX^ niV2 = MAX^ niVg = 0. 

(4) MAX^ niVi = {iyeNi\p = 2Kn}, 

MAX^ niVs = {iyeN2\p = 2Kn, cut 7^ 0}, 

MAX^niVg = 0. 

Proof. Apply Th. 16.11 Propositions 16. 3H6. 61 and Lemma [6^ D 



6.6 Upper bound on cut time 

The cut time for an extremal trajectory g^ is defined as follows: 

tcut — sup{ii > I gs is optimal for s G [0, ti]}. 

For normal extremal trajectories qs = Exp(A, s), the cut time is a function of 
the initial covector: 

tcut ■■ C ^ [0,+oo]. 

Denote the first Maxwell time as 

if A^(A) = inf{i > I (A,i) e MAX}. 

A normal extremal trajectory cannot be optimal after a Maxwell point, thus 

icut(A) < if^^(A) V A e C. 

On the basis of this inequality and results of Subsec. 16.51 we derive an effective 
upper bound on cut time in the sub-Riemannian problem on SE(2). To this end 
define the following function t : C -^ (0, +00]: 

A e Ci ^ t(A) = 2K{k), (6.18) 

AeCa ^ t(A) =2A:pJ(fc), (6.19) 

A e C3 ^ t(A) = +00, (6.20) 

A e C4 ^ t(A) = TT, (6.21) 

A e C5 ^ t(A) = +00. (6.22) 
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Theorem 6.3. Let XeC. We have 

icut(A) < t(A) (6.23) 

in the following cases: 

(1) XeC\C2, 

(2) A e C2 and snr 7^ 0. 

Proof. If A e Ci, then {X,4:K{k)) = (A,t(A)) G MAX^ by item (4) of Th. K^ 
thus 

<cut(A) < if ^^(A) < t(A). (6.24) 

If (A,i) G iV2 and p = p}(fc), snr ^ 0, then (A,i) e MAX^ by item (2) of 
Th. [13 and the chain ([eTM]) foUows. 

If A € C4, then the trajectories Exp(A,i) — (0, 0,— sit) and Exp(A*,t) = 
(0,0, Sit) intersect one another at the instant t — n, thus (A,7r) = (A,t(A)) G 
MAX'^, and the chain ((Oi)l follows as well. D 

6.7 Limit points of Maxwell set 

Here we fill the gap appearing in item (2) of Th. 16. 31 via the theory of conjugate 
points. 

A normal extremal trajectory (geodesic) qt is called strictly normal if it is 
a projection of a normal extremal At, but is not a projection of an abnormal 
extremal. In the sub-Riemannian problem on SE(2) all geodesies are strictly 
normal. 

A point qt of a strictly normal geodesic Qs — Exp(A, s), s e [0,i], is called 
conjugate to the point q^ along the geodesic qs'xiv — (A, t) is a critical point of 
the exponential mapping. 

It is known that a strictly normal geodesic cannot be optimal after a conju- 
gate point [3]. At the first conjugate point a geodesic loses its local optimality. 
Below we find conjugate points on geodesies with A € C2 not containing Maxwell 
points. These conjugate points are limits of pairs of the corresponding Maxwell 
points, the corresponding theory was developed in [16]. 

Proposition 6.7 (Propos. 5.1 [16]). Let v^, v'^ € N , Vn ^ i^'n, Exp(z^„) = 
Exp(t'^), n e N. Lf the both sequences {vn}, {'^nl converge to a point v = (A,i), 
and the geodesic qs — Exp(A, s) is strictly normal, then its endpoint qt = Exp(i') 
is a conjugate point. 

It is convenient to introduce the following set, which we call the double 
closure of Maxwell set: 

CMAX ={9eN\3{vr,^ (A„,t„)}, {v',^ = (A;,t„)} C N : 

Vn ^ J^^, Exp(z^„) = Exp(i/^), 71 e N, lim Vn = lim i/.'^ ^ i^\ ■ 
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It is obvious that Vn e MAX, thus CMAX c cl(MAX). 

Proposition 16.71 claims that ii v = (A, t) E CMAX and the geodesic Qs = 
Exp(A, s) is strictly normal, then its endpoint qt is a conjugate point. 

Proposition 6.8. Let v = {X,t) e N2 be such that p — pl{k), snr =0. Then 
the point qt — Exp(i') is conjugate, thus t > tcuti^)- 

Proof. Consider the points i^^ — {pl{k),T± 1/n, fc) e A'^2- Then i^+ 7^ i^^ and 
lim„_,oo !^^ = u. Formulas (|6.7p - (|6.1ip imply that Exp{u^) ~ Exp(i^+). Thus 
v e CMAX, and the statement follows from Propos. [6.7l D 

6.8 The final bound of the cut time 
Theorem 6.4. There holds the bound 

tcutW < t(A) VA e C. (6.25) 

Proof. Apply Th. [63] and Propos. EH D 

The function t(A) deserves to be studied in some detail. One can see from 
its definition (|6.18p - (|6.22p that the function t depends only on the elliptic co- 
ordinate k, i.e., only on the energy E (|4.ip of pendulum p.8p . but not on its 
phase ip. Thus we have a function 

t:^i-^t(^), t : [-l,+oo) ^ (0,+c»]. 

Proposition 6.9. (1) The functiont{E) is smooth for E E [—1,1)U{1, +00). 

(2) lim t(^) = tt; lim t(E) = +00; t - 2^-kNE ^ \ ^ Q as E ^ +00. 

Proof. (1) follows from smoothness of the functions K(}i) and p\(k) for k G 
[0,1). 

(2) follows from the limits lim K(k) = tt/2, lim K(k) ~ lim p](k) = 

+00, limp}(fc) = 27r. D 

fc— >o 

A plot of the function t{E) is given at Fig. [9l 

In our forthcoming work [19] we show that the inequality (|6.25p is in fact 
an equality, i.e., icut(A) = t(A) for X £ C. 
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Figure 9: Plot of the function E h^ t(_E) 
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